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Abstract. 

The trace of an arbitrary product of quantum operators with the density operator 
is rendered as a multiple phase space integral of the product of their Weyl symbols 
with the Wigner function. Interspersing the factors with various evolution operators, 
one obtains an evolving correlation. The kernel for the matching multiple integral that 
evolves within the Weyl representation is identified with the trace of a single compound 
unitary operator. Its evaluation within a semiclassical approximation then becomes a 
sum over the periodic trajectories of the corresponding classical compound canonical 
transformation. 

The search for periodic trajectories can be bypassed by an exactly equivalent initial 
value scheme, which involves a change of integration variable and a reduced compound 
unitary operator. Restriction of all the operators to observables with smooth non- 
oscillatory Weyl symbols reduces the evolving correlation to a single phase space 
integral. If each observable undergoes independent Heisenberg evolution, the overall 
correlation evolves classically. Otherwise, the kernel acquires a nonclassical phase 
factor, though it still depends on a purely classical compound trajectory: e.g. the fase 
for a double return of the quantum Loschmidt echo does not coincide with twice the 
phase for a single echo. 
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1. Introduction 

Notwithstanding the aptness of semiclassical (SC) approximations for uncovering 
classical structures underlying quantum evolution, their use for constructing ready 
algorithms to deal with increasingly complex experiments is yet to be established. The 
practical advantage of integrating Hamilton’s ordinary differential equations, rather than 
dealing with Schrodinger’s partial differential equation, is counterbalanced by the need 
to search for trajectories that are only indirectly specihed by boundary conditions, 
instead of arising directly from their initial conditions. This difficulty has led to 
the development of initial value methods (or inital value representations, IVR) that 
substitute the, so called, root search by an integration over families of initial valued 
trajectories p |2l [31 n El E]. Thus, IVR’s have been seen as a workable alternative, in 
spite of considerable criticism [7]. One of their main achievements is the evaluation of 
correlations for quantum operators, iiAB{t) [8], even though, if A is chosen as a density 
operator, this reduces to a single evolving expectation value {B{t)). Here we establish 
general SC approximations for multiple correlations of observables evolved by various 
unitary operators: (UBVC...). 

In a recent paper [9], henceforth labled I, the IVR approach to SC approximations 
was realized entirely within the Weyl representation, that represents the operator B 
by the phase space function R(x) = B{p,q), its Weyl symbol, or its Fourier conjugate, 
B{^), its chord symbol, another complete representation [10]. A remarkable feature 
of these particular representations (including the Wigner function [TT], in the case of 
the density operator) is that they are based, respectively, on reflection operators and 
translation operators [T2l fTS] ITO] . These are themselves unitary, so one can combine 
them with the evolution operators which act on each observable into a single composite 
entity. Then the expectation of an evolving observable was cast as a phase space integral 
over the Wigner function multiplied by the nearly classical function that represents the 
observable. 

This procedure is here generalized to the correlations of an arbitrary number, u, 
of observables undergoing general unitary evolutions. The evolving correlation depends 
on a single family of compound unitary operators, labeled by zz continuous parameters. 
For each parameter, the required trace of this compound operator is then obtained from 
the compound periodic orbits in the corresponding classical evolution, according to the 
standard SC procedure nms]. 

Even though this is an important step, the identihcation of the appropriate 
compound unitary operator does not free us from a search for orbits. It is true 
that continuous families of periodic orbits, within continuous families of canonical 
transformations, can be followed by a generalized Newton’s method, as the parameters 
are varied in small steps [16], but this is still a formidable task. Furthermore, such a 
reliance on continuity is at odds with the use of efficient Monte Carlo methods at the 
next stage, where one integrates over the parameters. It is then fortunate that the IVR 
approach can be extended to the general evaluation of correlations, by simply freeing 
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one of the parameters: The corresponding segment of the periodic orbit is then removed, 
so that one then deals with a reduced compound trajectory. This is still composed of 
evolutions intercalated by reflections, but now the trajectory is determined by its initial 
value. There always exists an extra reflection which closes such an open orbit, so that 
its reflection centre can be chosen as the extra free parameter. 

Just as in I, the IVR algorithm avoids caustic singularities, transforming them 
into nodal lines (or surfaces) of the compound propagator. There remains an overall 
ambiguity of sign to be determined as such a line is crossed, but the general procedure 
presented in na, henceforth labled II, can be immediately generalized for correlations. 

In the special case where all the evolution operators are metaplectic [^, the 
semiclassical theory is exact, including its IVR version. This provides scope for 
simple applications that illustrate the general features of the method, without gripping 
with the difficulties of a full SC calculation, as presented in I. Perhaps, the greatest 
simplihcation concerns caustics (or nodal lines in IVR): The important point is that, 
even though families of metaplectic propagators do cross caustics in any representation, 
as a parameter is varied, in their case, the hnal integral for the correlation has no risk 
of being divided into regions with different signs that need to be determined, as was 
discussed in II. 

The present paper follows closely on the track of I: The same notations are adopted 
and we incorporate here many relevant features. For instance, descriptions of any of the 
observables to be averaged may be supplied by the translation operators underlying the 
chord representation, instead of the reflections that belong to the Weyl representation, 
so that here we just focus on the latter. Again, we shall not develop explicitly the 
alternative of picking a pair of trajectories (forming a Final Value Representation, 
FVR) the possible advantages being discussed in I. These alternatives shall remain 
implicit so as to emphasize the purely original features of the present work. We shall 
also rely on the discussion of sign ambiguities associated with crossings caustics in II, 
since they are readily incorporated into the wider setting of evolving correlations. The 
great simplihcation here is one of scale: By focusing on mechanical observables for 
which the Weyl symbol coincides with the corresponding smooth classical phase space 
function except for corrections that are of hrst order in fi, one reduces the expression 
for the evolving correlation to a single phase space integral, irrespective of the number 
of observables. 

The following section presents the general construction of the appropriate 
compound unitary operator as the kernel for the correlation of evolving operators. 
Section 3 then interprets its SC approximation in terms of a compound canonical 
transformation dehned by a sequence of trajectory segments and derives its trace from 
the periodic orbits. The alternative IVR scenario is then developed in section 4, 
whereas section 5 presents the simplihcations inherent to the propagation of mechanical 
observables. All formulae are valid for an even number of observables. Modihcations 

§ Unitary operators corresponding to classical symplectic transformations, that is, linear canonical 
transformations, e.g. those driven by harmonic oscillators [UllIllEolEIlEaESlEllES]. 
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2. Compound unitary operators 


The outcome of a standard repeated experiment on a quantum system is expressed as 
an average over an (observable) operator, which may correspond to a standard classical 
variable, such as position, a projector, or a POVM. One can also measure correlations 
between such observables that have undergone different evolutions. In the simplest case, 
these may concern the same operator traversing coherently the alternative paths of an 
interferometer, or just measured at different times, as in the correlations of Leggett-Garg 
|26] . Then, so that the correlation is real, one evaluates some suitable symmetrization 
of 


C — {Aj^ity)... A 2 {t 2 )Aiifi)) — tr Ayity) ... A 2 {t 2 )Ai{ti) p, (2-1) 

where each of the Hermitian operators Aj{tj) undergoes a Heisenberg evolution driven 
by some unitary operator, Vj, that is 

Ai(ti) = dd As Vj{t,). (2.2) 


If one dehnes the intermediate steps as 

^i+i = (2-3) 

(with Ui = Vi and = 14 (G)^), the general form is obtained. 


C — tr Uu+iAi, Uy ... A 2 U 2 Ai Ui p, (2-4) 

in terms of the original observables Aj. No longer is one limited to symmetric Heisenberg 
evolutions, being that each observable can now be sandwitched between arbitrary 
unitary operators Uj. Hence, fl2.4D also includes evolutions such as the hdelity, i.e. 
the quantum Loschmidt echo [271 E]. i An example of direct application of such 
a time evolved correlation arises in the theory for time-resolved electronic spectra, 
depending on the evolution of two pairs of transition dipole operators. The Franck- 
Condon approximation then leads to an expression for the spectrum in terms of the 
hdelity, which was obtained in [28] using IVR. The present theory supplies in principle 
the full correlation witout any supplementary approximation. 

Following the same notation as in I and II, the centre symbol or Weyl symbol of 
operator A is 


A(x) = 2^tr 


Rx A 


(2.5) 


where N is the number of degrees of freedom, while the unitary operator, 
corresponding to the (classical) rehection through the phase space point x [T2[ [131110] , 
plays a fundamental role throughout. In other words, 24(x) is the Weyl representation 
of A. Alternatively, the chord symbol of the operator A can be dehned as 


A(|) = tr 


T t A 


( 2 . 6 ) 


II This is the case of a single observable Ai = I. 
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where is the Heisenberg operator corresponding to a phase space translation by the 
vector The chord symbol and the Weyl symbol are related by Fourier transformation. 
The advantage of both these representations is that their families of basis operators, 
{/?x} and {T^}, belong to the group of unitary operators. An arbitrary operator is then 
expressed as a superposition of unitary operators 

dx ... - 1 


A = 2 


N 


A(x) = 




. {27rh)^ ^ {27rh)^ 

though it is convenient to keep the special notation for the density operator, 


p = 2 


N 


(2.7) 


( 2 , 8 ) 


dx I-F(x) = jdi x(i) 

in terms of the Wigner function [11], hF(x), and the chord function [10], x(^). 

In the case of the Weyl-Wigner representation, one can now insert these in the 
expression for the evolving correlation 

2N . 


C = 


dxi,...(ix 2 dxidxo A,,(x,,)...A 2 (x 2 ) Ai(xi) IF(xo) tr U{xo, xi, X 2 , ...,x,,}, (2.9) 


where the family of compound unitary operators is defined as 

U{xq, Xi, X2,..., Xjy} = Uy... Rx^Ui Rxq' (^'10) 

Before any evolution takes place, that is, when all of the tj = 0 in fl2.ip . the 
compound operator is just a product of reflections, because all the Uj = I, the identity 
operator: 

Uo{Xo,Xi,X2, ...,Xj,} = .Rx. ••• Rxi Rxo- (2-11) 

The simplest case is when u, the number of reflections is odd, i.e. an even number 
of observables. Then the product is also a reflection and we can identify tr Uq with 
the kernel for the Weyl representation of a product of Weyl operators, each of which is 
specihed by its Weyl symbol [TO] : 


tr Uo{xo, xi, X 2 ,..., X,,} = 2 


-uN 


exp 


-Aj,+i(xo,Xi,X2, ...,x^ 
M 


( 2 , 12 ) 


where Ay^i is the symplectic area of the polygon whose sides are centred on 
{xq, Xi, X 2 , ..., Xjy} as drawn in Fig 1; this is a bilinear function of each pair of variables, 
which arises in the general product formula for the Weyl representation of the product 
of an even number of operators, Ay...Ai, in [10], that is. 


{A^...Ai}(xo) = J A^(x^)...Ai(xi) exp 


-A^+i(xo,xi, ...,x^) 
n 


(2.13) 


{nny^ 

so that one retrieves the simple expression for the initial correlation: 

Co = Idxo {A,...Ai}(xo) IF(xo). (2.14) 

In the following section, we show how this polygonal scenario is extended to 
evolving correlations within the SC approximation. Evidently, one can express evolving 
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Figure 1. The initial kernel of the integral for the correlation of an even number 
of observables is the complex exponential of Ai,_|_i(xo, xi, X 2 ,Xi,), the symplectic 
area of the unique polygon with sides centred on {xq, xi, X2,x,y}. For, v = 2, 
this is a triangle (left), whereas the plane projection of higher polygons may have 
self-intersections (right). 


correlations equally well in terms of compound operators with translations in place of 
reflections, 



(2.15) 


or it may be more convenient to keep some reflections and some translations. The issue 
appears in the context of a single evolving expectation and it is already discussed in I. 
Indeed, the product of any even number of reflections, Uo{xo, xi, X 2 ,..., Xj,}, results in a 
translation operator, rather than a reflection, so that it is advantageous to also use the 
chord-translation basis for the density operator. For this reason, the results that must 
be adapted for an odd number of observables, by resorting to the Fourier transform of 
the Wigner-Weyl representation, will be remitted to the appendix. 

3. SC approximation for the compound Weyl propagator 

The key ingredient for the SC approximation of the compound propagator and its trace 
is the general SC Weyl propagator corresponding to an arbitrary classical canonical 
transformation that is generated by a Hamiltonian, H(x.), acting during a time, t. In 
the simplest case, this is simply [291 HU] 



(3.1) 


The geometric part of the centre or Weyl action, S'(x), is just the symplectic area 
between the trajectory and the chord, ^ — x", joining its endpoints. From this, 

one subtracts —Et, where E is the energy of the trajectory. The Maslov index, a, is zero 
in a neighbourhood of the identity operator, in which case there is indeed only a single 
classical trajectory centred on the point {x = (p,q)} of classical phase space 
Otherwise, there may be multiple solutions to the variational problem that identihes 
trajectories with a given centre, x, so the actions may have many branches and these 
branches meet along caustics where the semidassical amplitude diverges. (See II for 
the phase a and further details.) 
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The centre action specifies the classical canonical transformation, x x+, 

corresponding to U indirectly through [10] 

« = X+ = X+|. x-=x-|. (3.2) 

The linear approximation of this transformation near the x-centred trajectory is defined 
by the symplectic matrix M. This has the Cayley parametrisation: 

M = [1 +JB]-^[I-JB], (3.3) 

in terms of the symmetric matrix B, which is just the Hessian matrix of «S'(x). 

A notable exception for this SC form of the Weyl propagator is precisely that of a 
reflection operator, i?x- Indeed, its Weyl symbol is simply, 

i?x(x') = 2-^5(x'-x). (3.4) 


It is the chord representation of this operator that has the standard semidassical form. 

To construct the SC approximation of the compound propagator, U{xo, xi, X 2 ,..., Xj,}, 
we assume that each of the Weyl propagators, ie. the Weyl symbols, Cj (x), for Uj can 
be expressed in the form fl3.ip . Then the key point is that the compound unitary op¬ 
erator, U{xo, xi, X 2 ,..., Xj,}, for any choice of parameters, has its own Weyl symbol, 
U{xo, xi, X 2 ,..., Xj,}(x) or U(x) for short. This compound propagator is determined 
from the factor propagators, hj(x') according to the product formula (I2.13p : 


U{xo,Xi,X2, ...,Xj,}(x) = 


(ix(,,;^dx"...(ixjx; 


(7rfi)2C+i)^ 


c.+i(x;+i) i?x.(x'')...f/i(x;) i?xo(x'') 


exp 


A ' / ' 

2^+3(Xo,Xi,..., 




(3.5) 


The integrals over the arguments of the reflections merely fix the respective centres, 
x" = Xj, whereas all other integrals are evaluated semiclassically by stationary phase. 
But the deduction of (13.1 p in [10] was itself based on the same product formula, that 
is, from a product of infinitesimal small time propagators, so the result is again a Weyl 
propagator of the same form, but it is constructed from a compound trajectory built 
up from the sequences of partial trajectory segments. The relevant trajectory, which 
is centred on x, i.e. the argument of the Weyl propagator, is built up out n + 1 such 
segments joined hy n + 1 reflections. 

The appropriate trajectory segments that are generated by each centre generating 
function, Sji^x'j), as well the centres themselves, need to be chosen so that the full 
compound trajectory is continuous. This is achieved by imposing the requirement [TO] 
that the overall action in (13. ip is just 


S(x) — A21/+3 -l- i 5 'i(x'i) -|-... -|- S'i,+i(x',^+i). 


(3.6) 


Here A 21.+3 is the symplectic area of the dynamical polygon with a side centred on x, 
as well as u + 1 sides centred on the points Xj and u + 1 sides centred on x'^. The 
stationary conditions. 


as _ ^ dA2.+2 _ ds, _ . 

ax', ° dx', dx', 


(3.7) 
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define the variables x'j, so that the trajectory arcs ht precisely the sides of the dynamical 
polygon, as depicted in Fig. [2l which exemplifies the general layout for the case of two 
observables, i.e. n = 2. It is naturally satished by any compound trajectory that is 
followed through from its initial value. As for the compound monodromy matrix, M, 



Figure 2. Phase space structure corresponding to the compound Weyl propagator 
for the evolution of a pair of observables. The arguments of their Weyl symbol are the 
reflection centres xi and X2, while the reflection centre for the Wigner function is Xq. 
These reflections are joined by the trajectories of the intervening SC Weyl propagators 
with centres at x'l and x '2 so as to form the dynamical polygon. (These trajectories 
need not be generated by the same Hamiltonian) 


this is just the product of the sequence of monodromy matrices for each step. Indeed, 
since the matrix for a reflection is just —I, we have simply 

M = [-1] ■ Ml • [-1] • M 2 • ...[-I] ■ M,+i = (-l)"+iMi ■ M 2 • ...M,+i.(3.8) 

The kernel for the propagation of the correlation between observables, Aj, 
evaluated at multiple times in terms of the initial Weyl symbol, Aj(x), is simply 
tr U{xo, xi, X 2 ,..., Xy}. Just as for other representations, the SC limit of the trace 
of an unitary operator singles out the contributions of the periodic trajectories. In the 
case of the Weyl symbol, this is derived from 

‘^® = /( 2 ^U(x). (3.9) 

Now, the only explicit dependence on the argument, x of the compound propagator is 
in the symplectic area of the dynamical polygon, A 2 n +3 = x A ^ + Const, where $, is the 
side centred on x. Like the term Const, this depends only on the other centres [lO] : 


Thus, 


— = (x 1 
2 ^ ^ 


(Jx 


{27rH) 


N 


Xo) 


exp 


(X2 - XiJ 


(X ^+1 - Xj, 


-A 


2^+3(X,Xj,XjJ 


= exp 


:A2^+2(xj,x'j) 


(3.10) 


<5(0. (3.11) 


that is, the integral singles out those polygons where the open side has zero length. The 
compound action for the trace is the same as 


), but with A 21/+3 A 21 .+ 2 , which is 

the polygon for a periodic orbit. Such a compound periodic orbit for the trace is depicted 
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Figure 3. The SC expression for trace of the compound propagator is determined 
by the corresponding periodic trajectories: The open side of the dynamical polygon, 
A 7 in the previous figure collapses, so that x = x“ = x+. 


in Figini The corresponding monodromy matrix is just fld.Sp without the last factor, 


l-Il 


The stationary condition for the intermediate centres, xh, of the periodic orbit 


is just fl3.7p . so that this merely prescribes each side of the dynamical polygon to be 
the chord for the corresponding trajectory segment. However, in practice one need not 
search for the solution of each of these equations because they are automatically satished 
for any given periodic orbit of the compound canonical transformation, i.e. one then 
has all the centres and the sides of the (2i^)-sided polygon for a particular pinning of the 
reflection centres {xo,xi, ...,Xy}. Then one needs only to evaluate (13.911 by stationary 
phase, so that the contribution of the p’th periodic orbit to the trace is just 




tr Up 


det B det(I 

2N 


M)|V2 


exp 


Uo/ \ ti'^d. 

- S X + 


det(I - M)|V2 

where, following [10], one takes the determinant of (13.3p . 0 The parameters 

{xo,xi,X 2 , vary continuously in the ultimate integration for the correlations, so 

that the periodic orbits at each neighbouring parameter can be found by a generalized 
Newton’s method, such as in [T6] . 

The amplitude of the periodic orbit contribution to tr U{xo, xi, X 2 ,..., Xj,} becomes 
singular at those values of {xo,xi,X 2 , ...,Xy} that determine periodic orbit resonances: 
det (I — M) = 0. This problem is avoided by the generalization of the IVR theory of I 
in the following section. The IVR approach even dispenses with the search for periodic 
orbits in the evaluation of evolving correlations. 


4. Initial value representation 

Let us reinterpret the kernel of the evolving correlation (12.911 as 

tr U{xo,xi, ...,x^} = tr Rx. Vj.... RxiVi Rxo = U'{xi,..., x^}(xo), (4.1) 


^ The stationary phase evaluation adds a further phase of 7 r /4 times the signature of B. This will not 
be needed in the IVR theory in the next section, so that it is here just included in the Maslov index a'. 
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that is, according to fl2.5p . the Weyl symbol for the reduced compound operator. 

U'{xi, = U^+i U^... R^^Ui. (4.2) 

In its terms one generalizes the Weyl representation to the evolved product operator in 
(I2.13P as 

2^ f 

{Aj....^i}'(xo) = J dxj....dxi U'{xi,Xj.}(xo), (4.3) 

so that the correlation is simply 

C = I d^o f^(xo) {2l,...2li}'(xo). (4.4) 


Thus, U' has one less reflection than U, but clearly its SC approximation has 
the same form as fl3.ip . with the corresponding classical polygonal path given by 
Fig 121 rather than FigJ2] . In other words, FiglHl is now reinterpreted as an open 
polygonal line, with its begining and end points, xq^, centred on xq, as shown in 
FigJU It is important to note that, even though the full canonical transformation, 
Xq” xq^, can be decomposed into several partial canonical transformations, the 
quantum unitary transformation U' corresponds to this single reduced compound 
canonical transformation. Furthermore, each branch of its centre generating function , 
S'(xo), in the SC approximation to the Weyl propagator, U'(xo), is constructed from 
those compound trajectories that satisfy xo“ + xq^ = 2xo. 



Figure 4. The same polygonal trajectory as in the previous figure is reinterpreted as 
an open trajectory of the reduced compound canonical transformation, from which the 
initial reflection at Xq is removed. This point is now the argument of the Weyl symbol 
for the reduced compound propagator. 


The initial value representation (IVR) now results from the exchange of the 
integration variable from the trajectory midpoint, xq, that is the argument of the Wigner 
function in the evolving correlation (12.hh . to the initial point, xo“. This generalization 
of the procedure in I relies on the simple form of the Jacobian of the transformation: 


detih^ = detil±^' 


d^o 

Thus, it brings the SC approximation for the evolving correlation to the form 

1 


(4.5) 


C 


(ttH) 


uN 


y dxi,...(ixi(Jxo Aj,(xj,)... Ai(xi) VF(xo(xo ))| det(I + M')|F2 
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exp 


;(S'(xo(xo )) + nna) 


where M' is the monodromy matrix for the linearization of xq 
neighbourhood of xo“. It can be decomposed as the product 


(4.6) 


H->■ Xq’*' in the 


M' = Ml • [-1] ■ M2 • ...[-I] ■ M,+i, (4.7) 

a reduced version of fl3.8p , but now the sequence of factor monodromy matrices is directly 
determined by the initial value xq”. Likewise, each of the variables x'j will be just the 
centre of of the respective side of the polygonal path starting at xq”. The crucial 
point is that the Jacobian for the switch to the new integration variable, xq~, kills off 
the caustic singularities in the SC kernel for the evolving correlation. Thus, in a single 
step, without increasing the number of integrations, one does away with the need to 
search for trajectories while erasing all caustics! 

The generating function, S'(xo), for the canonical transformation is now defined as 

S'(xo) = A'2,.+2(xo, X 1 ,..., x^, x^+i) + 5'i(x'i) + ... + 5'i.+i(x'^+i), (4.8) 


where the requirement 

<9S' ^ dA'2.+2 dS, ,, 

defining the variables x'j is automatically satisfied by any compound trajectory that is 
followed through from its initial value. Indeed, as discussed in [10], the symplectic area 
of the reduced dynamical polygon satishes 

2A'2i/+ 2 = ^\A^]^ + (^\+^]^)A|'2 + --- + (^\+^i+^^2+-'-+^^i/)/\^^iy+i) (4T0) 



where each chord refers to the appropriate reflection (unprimed) or partial evolution 
(primed). The discussion of the Maslov phase a in II applies direcly to fl4.6p . 

In a full SC calculation in which each trajectory segment needs to be integrated 
numerically, the numerical error will build up along the sequence of segments. It may 
then be preferable to start somewhere in the middle of the sequence, from where the 
sequence is taken partly backwards and partly forwards. This is a direct generalization 
of the Final Value Representation (FVR) for the evolved average of a single observable 
that was presented in I, even though the forward and backward paths only have the 
same number of segments if v is odd. In any case, the Jacobian for the exchange of the 
initial value for an intermediate value, xq” ^ Xj“, a canonical transformation, is just 
unity. 


5. Evolving mechanical observables 

So far no account has been taken of features that distinguish observables from other 
operators. Within the Weyl representation, mechanical observables are represented by 
real smooth functions of the phase space variables. Indeed, the Weyl representation 
of the operator function, A(x), is just the classical phase space function, A(x), plus 
corrections of order ft, which depend on the chosen symmetrization of products of p 
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and q. Even more to the point, the Weyl representation of a prodnct of mechanical 
observables, Ai^...Ai is just {Ai^...Ai}{'x.) = 74,^(x)...74i(x), up to first order terms in H, 
which again depend upon ordering. 

It is important to understand how this simplification arises, starting from fl2.13p . 
So one adapts the discussion concerning equations (Id.llh and (14.91) : In the absence of 
any other phase term beyond the symplectic polygonal area, stationary phase evaluation 
of the multiple integral in fl2.13p for all the variables {xj,...xi} collapses each side of the 
polygon in Fig.l, = 0 ... = 0, and hence the polygon itself, with Xj, = ... = xi = xq. 

The fi-dependent corrections can be calculated via a generalization HDl of the familiar 
Groenewold product formula inni, but they will be only of second order in h for 
symmetrizations that guarantee a Hermitian product. 

This simple classicality of the Weyl representation of a product of mechanical 
observables, which is shared by the initial correlation, 

Co ^ y dxo A^(xo)...Ai(xo) W(xo), (5.1) 

may be destroyed as the observables evolve. Even so, the absence of high period 
oscillations in each of the Weyl symbols, Aj{x.), still allows for stationary phase 
evaluation of fl4.3p . analogous to that of fl2.13p . Indeed, the expression fl4.8p for 
the reduced action, S"(xo), can now be reinterpreted as the composition of 2;/ + 1 
transformations, but with zero action for the unprimed variables. El Then the stationary 
condition (14.91) for each of these centres collapses the corresponding side of the (reduced) 
dynamical polygon, that is ^j(xj) = 0, as depicted in Fig 5. Nonetheless, the evolution 



Figure 5. Stationary phase integration over the unprimed reflection centres collapses 
the corresponding sides of the dynamical {2u + l)-sided polygon: These corners of the 
resulting [v + l)-sided polygon now depend on the initial value, Xo“. 

sides, ^h(x') are no longer zero, so there is generally a non-zero phase in the integrand 
that is constructed from the remaining {p + l)-sided polygon: Each unprimed variable, 
Xj is now a free corner, depending on the initial value xq”, instead of being the fixed 
centre of a side. In this way there results an enormous simplification of the correlation 

+ Note the subtle difference with respect to the multiple integral (13.51) : There, one integrates over the 
primed variables, arguments for the Weyl symbols for Uj. Here, we integrate first over the unprimed 
variables, which represent the mechanical observables, Aj. 
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formula 04.61) : 



/ 


dxo ))... Ai(xi(xo )) ll/"(xo(xo ))| det(I + 



h 


(5.2) 


Thus one no longer deals with a full family of compound canonical transformations. 
Instead of this, the polygonal trajectories are built for each initial value within a single 
canonical transformation without the intermediate reflections. 

In the limit of short times, u sides of the remaining [n + l)-sided polygon shrink 
to zero. In this limit the overall phase is zero, so that one retrieves (15.Ij) if v is even. 
This restriction on n follows from the expression fl3.8p for the monodromy matrix in the 
amplitude of the compound propagator: As each of the matrices Mj I, one obtains 
det(I + M') —)■ 2^^ if n is even, or zero if n is odd (and hence a caustic). In the latter 
case, the Appendix obtains the correlation as a single integral of the Weyl symbols for 
the observables weighed by the chord function instead of the Wigner function. 

The case of multiple Heisenberg evolution fl2.2p also collapses the SC phase, but for 
all time! The easiest way to see this is to propagate directly the Weyl representation of 
each observable, A{t): 



so that the phase space point representing A is placed at the corner of the triangle 


oposite the point x where the evolved observable is evaluated. In HU it is shown that 
A(x, t) = A(x'(f,x)), the classically evolved observable, for a metaplectic evolution, i.e., 
for a quadratic driving Hamiltonian. This is also the correct semiclassical approximation 
for a mechanical observable, resulting from stationary phase integration: There is only 
a single trajectory traversed both forwards and backwards from the initial point, x, and 
Xi = X 2 is its midpoint, so that there is complete phase cancellation. In the case of 
the full correlation, the curved polygon in Fig 5 for the correlation collapses into a thin 
legged u-spider as shown in Fig 6. Then one can merely obtain the evolved correlation 
from fl5.ip with the classically evolved observables: 



(5.4) 


It is important to note that the collapse of the dynamical polygon into a spider concerns 
exclusively the pairing of classical trajectories corresponding to the unitary operators 
within the reduced compound operator. Thus there is no restriction on the Wigner 
function, that is, the evolution is purely classical even for a highly oscillatory quantum 
Wigner function. 

6. Discussion 


The elaborate theory underlying SC approximations for the multiple evolution of 
correlations between arbitrary numbers of quantum observables leads to deceptively 
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^2 

Figure 6. The dynamical polygon collapses for multiple Heisenberg evolution 
of mechanical observables, with each point Xj(xo,t) placed at the tip of a classical 
trajectory that is retraced with negative time. The spider shape arises when each 
observable is propagated by a different Hamiltonian. There is no phase factor in this 
essentially classical evolution for the correlation between mechanical observables. Any 
difference between forward and backward pair of Hamiltonians fattens the legs and 
these separate, even if the same pair propagates all the observables. 


simple results. The approximate semidassical scenario presented here becomes exact 
in the limit where all the evolution operators are metaplectic, that is, generated by 
quadratic Hamiltonians. 

Correlations are independent from the choice of representation that is employed in 
their calculation, but the reliance on Weyl symbols and the Wigner function for the 
present theory leads to a rich phase space structure. The evolution kernel is identihed 
with a compound propagator, which is constructed by a dynamical polygon whose sides 
are orbit segments that are combined into a compound classical trajectory and thence 
to a quantum phase. Depending on the choice of symmetrization of the observables, an 
appropriate trigonometric function of the polygonal area will dephase the correlation 
integral. 

It is shown in the Appendix that the need to distinguish whether the number 
of observables is even or odd dissolves for short times and for the restricted class of 
mechanical observables, such that the Weyl representation of their product coincides 
with the corresponding smooth classical phase space function, except for corrections 
that are of hrst order in ft. The general rule is that independent Heisenberg evolution 
for each observable leads to fl5.4p a classical expression of the evolving correlation, even 
if the Wigner function employed in their average has quantum oscillations that are 
separately though concurrently sampled by each observable. 

In contrast, if each observable does not follow its own Heisenberg evolution so that 
the intermediate evolution operators are not given by fl2.3h . a phase factor will grow 
in time within the single phase space integral for the correlation. This is generalizes 
the result for the hdelity (or the quantum Loschmidt echo) presented in I, a special 
case within the present framework, that of a single observable, the identity operator, 
undergoing different forward and back evolutions. In the language of section 5, the 
dephasing grows with the symplectic area of a single slightly fattened spider leg, i.e. a 
curvilinear triangle. What about a repeated echo: On returning, one evolves again with 
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the same pair of forward and back Hamiltonians? Then the initial value for the second 
traversal has changed, so that a new spider leg is drawn which is only initially close to 
the hrst leg, even though it is generated by the same pair of Hamiltonians (specially if 
they are chaotic). The correct symplectic area that determines the dephasing is then 
that of the full two legged spider, a curvilinear pentagon, instead of twice the area of the 
first triangular leg. For small times, the difference may be small, but the denominator 
for the phase is Planck’s constant... 

Whereas the simple dephasing representation of Vanicek inu atributes the dephasing 
factor for the hdelity to a single classical trajectory, this was shown in [32] to result from 
a first order classical perturbation theory of an action for a pair of trajectories. A further 
generalization in I related the evolution of the expectation of a single observable to a 
pair of trajectories surrounding a translation or a reflection. Now we hnd that the only 
price to pay for having more mechanical observables in a correlation is to add more 
segments to the corresponding compound classical trajectory. In all cases, the relevant 
trajectory is completely specihed by its initial value, i. e. the integration variable in 
the average. Furthermore the general analysis of Maslov phases in II, that is valid for 
all cases, guarantees that initially, they are absent and it is only after a first caustic is 
crossed that extra phases need to be taken into account. 

One should note that relaxing the restriction to mechanical observables does 
not necessarily complicate matters. Observables may well be projectors, so that 
the correlations become correlations between measurements. For instance, one may 
measure the momentum (see [8] for examples), i.e. the projector. Bp, rendered in the 
Weyl representation by the singular distribution, i?p(x) = 5{p — P), which actually 
simplihes the multiple integral for a correlation. More generally, positive operator 
valued measures (POVM, see [18]) are also represented by phase space functions in 
the Weyl representation. A specially interesting case to study in this setting is that of a 
parity projection onto the subspace for either eigenvalue ±1 of the reflection operators. 
Reflections are observables as well as being unitary operators and their Weyl symbol is 
just a pointwise delta function. Their measurement was proposed in [33l [3l] and carried 
out experimentally in [35] . Thus, one may readily extend the present results beyond the 
restricted class of mechanical observables. 

Appendix A. Mixed centre-chord propagation kernel 

Before any evolution takes place, that is, when all of the tj = 0 in (I2.1[l . the compound 
unitary evolution operator is just the product of reflections fl2.1ip . When n, the number 
of reflections is even, i.e. an odd number of observables, the product is a translation, 
rather than a reflection and its trace is a delta function in the Weyl representation ng. 
Thus the zero time limit of tr U{xo,xi, = U'(xo) is singular. 

On the other hand, one may resort to the chord representation of the reduced 
compound unitary operator, 

tr Rx. Uu--- Rxi^i T_^ = U'{xi, ...,Xj,}(^o), 

So 


(A.l) 
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even though it is specihed by its intermediate reflection centres, so that the evolved 
product operator has its chord representation, 

j ^^(x^)...v4(xi) U'{xi,...,x^}(^o), (A.2) 

similarly to its Weyl representation fl2.13p . Thus, the general expression for the trace of 
a product in the chord representation [10] supplies the correlation as 

C = x(-«o), (A.3) 

whether or not the operators Aj have evolved. 

In effect this is a generalization of the correlation for a single observable, i.e. its 
expectation in I: The evolution was there atributed to the density operator, rather 
than to the single observable, but here the opposite choice is preferable, because the 
observables may evolve independently. One should note that, in the case of an even 
number of observables, it is the chord representation that is singular initially. So there 
is no way around the need for different treatments depending on the parity. 

To obtain the initial correlation for an odd number of observables, one needs the 
chord symbol for an odd number of reflection operators. Following the relations in ra. 
this is deduced to be 



Xl 




'XO 5 


(A.4) 


where A,^ is the area of the closed polygon with sides centred on {xi, ...,x,x} as shown 
in Fig.l, even though the polygonal path generated by u reflections from an arbitrary 
point is not closed. The segment that does close such a. {n + l)-polygon is necessarily 
centred on 

xo = Xj. - t^^(xi, ...,Xj._i), (A.5) 

where, in its turn, is the open side of the i^-polygon whose other sides are centred on 
{xi, ...,X^_i}. 0 It follows that 




(A,6) 


SO that the chord symbol for an arbitrary product of operators, each one specified by 
its Weyl symbol, is 




dxj,...(ixi 


Ai.(xj.)...Ai(xi) 




(A.7) 


Thus, according to flA.3p . the initial correlation in the case of an odd number of operators 
is 


Cn = 


(ixy...(ixi j ^ ^ X / \ -^Ay(xi,...,xi.) f d^o X 


J .. / ^ x(-«„) et 

/ "’‘‘'V(xo(xi, ...x^)). (A.8) 


* If = 3, xo is just the missing corner of the parallelogram with its other corners at {xi, X 2 , X 3 }: The 
inscribed polygon defined in uni- Inscribed polygons, for all odd v, satisfy special constraints. 
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If all the operators, Aj, are now assumed to be mechanical observables with nearly 
classical, non-oscillatory Weyl symbols, Aj(xj), then the stationary phase evaluation 
of the integral over Xj collapses the side of the polygon, which it centres, i.e. = 0. 
Performing {v — 1) stationary phase integrations sequentially, from Xj, to X 2 , the polygon 
looses its sides, while the argument of the Wigner function in the remaining integrals 
becomes dependent on fewer variables. Finally, the last remaining integral is just 

Co ^ y dxi A^(xi)...Ai(xi) IF(xi), (A.9) 

which coincides with (ED, inspite of the number of observables being odd. This may 
now be reinterpreted, such that the Weyl symbol for the product is approximately 
A,^(x)...Ai(x), so that the symbol for the evolved product is {Ay.AiYAA ~ 
ylj,(x)...Aj(x) and hence that (15.41) holds approximately for the independent Heisenberg 
evolution of all operators, whether even or odd. 
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